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Abstract 

Recently, the ATLAS and CMS collaborations have announced the discovery of a 
125 GeV particle, commensurable with the Higgs boson. We analyze the 2011 and 
2012 LHC and Tevatron Higgs data in the context of simplified new physics models, 
paying close attention to models which can enhance the diphoton rate and allow for 
a natural weak-scale theory. Combining the available LHC and Tevatron data in the 
h -)• ZZ* ->■ 41, h ->• WW* ->■ Mu, h ->■ 77, hjj ->■ jjjj and hV ->• bbV channels, 
we derive constraints on an effective low-energy theory of the Higgs boson. We map 
several simplified scenarios to the effective theory, capturing numerous new physics 
models such as supersymmetry, composite Higgs, dilaton. We further study models 
with extended Higgs sectors which can naturally enhance the diphoton rate. We find 
that the current Higgs data are consistent with the Standard Model Higgs boson and, 
consequently, the parameter space in all models which go beyond the Standard Model 
is highly constrained. 



1 Introduction 



ATLAS and CMS [JJ [2] have just presented an update of the Higgs searches, independently 
combining about 5 fb _1 of data collected at = 7 TeV and more than 5 fb" 1 at = 8 TeV. 
Both experiments observe a distinct excess in the diphoton invariant mass spectrum near 
125 GeV with a local significance of 4.5a and 4.1a, respectively. Moreover, an excess of 4- 
lepton events with my ~ 125 GeV, which can be interpreted as a signal of the h — > ZZ* — > 41 
decay, is observed by both experiments with the significance of 3.4a and 3.2a, respectively. 
Combining all available channels, the significance of the signal is around 5.0a for both 
ATLAS and CMS. Given these data, the existence of a new resonance near 125 GeV is now 
established beyond reasonable doubt. The remaining open question is the precise nature of 
that resonance. Is it a Higgs boson? Is it the Standard Model Higgs boson? If not, what sort 
of new physics is being favored or disfavored by the Higgs data? 

To address these questions, in this paper we combine the available ATLAS [3J HJ |5], E] , 
CMS [TJ, [HI El ID] , and Tevatron [H] data in several channels that are currently most sensitive 
to a 125 GeV Higgs. We interpret the data as constraints on an effective theory describing 
general interactions of a light Higgs boson with matter [T21 [T3"| H3] ■ in this approach, a small 
number of couplings captures the leading-order Higgs interactions relevant for production and 
decay processes probed by the LHC and Tevatron. The Higgs of the Standard Model (SM) 
corresponds to one point in the parameter space spanned by the effective theory couplings. 
Given the event rates measured by experiments and the corresponding errors (assumed to 
be Gaussian), we can construct the likelihood functions in the parameter space. This allows 
us to determine whether the SM Higgs is consistent with the existing data, and quantify 
the possible departures from the SM in a general framework with a minimal number of 
theoretical assumptions. 

We also confront the Higgs data with a number of simplified models beyond the SM (for 
earlier such studies based on 2011 data, see e.g. [15] and reference therein). New degrees 
of freedom coupled to the Higgs and carrying color and/or electric charge may affect the 
effective couplings of the Higgs to gluons and photons, while mixing of the Higgs with 
beyond Standard Model (BSM) scalars may affect the Higgs couplings to the W, Z and the 
SM fermions. In this context, we discuss which couplings and mixing angles are allowed 
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by the data. Unsurprisingly, preferred models feature an enhanced rate in the diphoton 
channel, as indicated by the data. We further put a special focus on whether the allowed 
regions are consistent with natural theories, where the new degree of freedom cancel the 
quantum corrections to the Higgs mass induced by the SM top and electroweak bosons. 

The paper is organized as follows. In the next section we define the Higgs low energy 
effective action and identify the relevant parameters that are being constrained by the present 
data. We calculate the contributions to the action from integrating out new physics particles 
and derive the mapping of the action to the corresponding rates which are measured by 
the experiments. In Section [3] we discuss the LHC and Tevatron Higgs data and show 
the resulting constraints on the parameters of the Higgs effective action. In Section [4] we 
then study the constraints on simplified models. In doing so, we capture several motivated 
scenarios such as little Higgs and supersymmetric models, W, and dilaton. Section [5] is 
devoted to models with extended Higgs sectors. We show that simple extensions such as a 
doublet-singlet or doublet-triplet Higgs sectors, allow for an enhancement in the diphoton 
rate in agreement with all other constraints. We conclude in Section [6j 

2 Formalism 

We first lay out in some detail our effective theory approach. We describe interactions 
of the Higgs boson with matter using an effective Lagrangian where a small number of 
leading order operators capture the salient features of Higgs phenomenology. Using the 
effective Lagrangian we derive the relevant production and decay rates as a function of the 
effective theory couplings. With these relations at hand, one can then construct the coupling- 
dependent likelihood function for a set of measurements, allowing for bounds to be placed 
on these couplings and the best-fit regions to be identified. We closely follow Ref. [12J; for 
similar approach to analyzing the 2011 Higgs data, see [131 [H]. 
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2.1 Higgs Effective Action 

We introduce the effective Lagrangian defined at the scale of \i = m h ~ 125 GeV, 

1m\? „,,„._ , m\ m b U Ti m r u - m cu- 

*->eff = Cy ^W^W + C V hZnZn — Cb ll 00 — C T II TT — C c (X CC 

V M M V V V V 

QL ol 

+ c gT7T—hG a G a + ay—hA^A^ - c inv h\X • (2.1) 

This Lagrangian describes the interactions of a light Higgs scalar with matter, providing a 
very general and convenient framework for interpreting the current Higgs searches at the 
LHC and TevatronQ The couplings of the Higgs boson are allowed to take arbitrary val- 
ues, parametrized by To be even more general, we also allow for a coupling to weakly 
interacting stable particles \i leading to an invisible Higgs partial width [TU [TB] . This effec- 
tive approach harbors very few theoretical assumptions. One is that of custodial symmetry, 
cw = cz = cy so as to satisfy the experimental bounds on the T-parameter, with h assumed 
to be a singlet of custodial isospin. Another theoretical assumption is that the Higgs width 
is dominated by decays into up to 2 SM particles; more sophisticated BSM scenarios may 
predict cascade decays into multiple SM particles which would require a separate treatment. 
Finally, we assume that the Higgs is a positive-parity scalar; more generally, one could allow 
for pseudo-scalar interactions. 



The top quark has been integrated out in Eq. (2.1) and its effects are included in the 
effective dimension-5 Higgs couplings parameterized by c g and c 7 . However these 2 couplings 
can receive additional contributions from integrating out new physics, and therefore are also 
kept as free parameters. At the same order one could include the dimension-5 Higgs coupling 
to WW and ZZ, however their effects can be in most cases neglected in comparison with 
the contribution proportional to cy. We therefore omit them for simplicity. The Lagrangian 
should be extended by the dimension-5 coupling to Z7, once measurements in this channel 
become available [TT] . To describe the tt associated Higgs production process, which may 
be observable in the 14 TeV LHC run, one would not integrate out the top quark. 



1 A tacit assumption is that we are using the effective Lagrangian to study processes where the Higgs 
boson is dominantly produced near threshold. For exclusive processes requiring Higgs produced with a large 
boost, PT,h S> mh) the contribution of higher order operators may be quantitatively important. 
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2.2 SM and New Physics Contributions 



In the SM, the terms in the first line of Eq. (2.1) arise at tree-level: 



CV,SM — Q),SM — C T; sm — C Ci SM — 1- (2.2) 

The values of these couplings may be easily affected by new physics if the Higgs boson mixes 
with other scalars, or if the SM fermions and gauge bosons mix with new BSM particles. As 
we show below, one may easily construct models where these couplings are either enhanced 
or suppressed. 

We now discuss the possible contributions to the dimension-5 Higgs couplings in the 



effective Lagrangian (2.1). Consider a complex scalar S, a Dirac fermion /, and a charged 



vector (assumed colorless for simplicity) coupled to the Higgs as 

£ = _ c ?™Z h tfs - cp-hff + c^fcpk. (2.3) 



V V V 

Integrating these particles, the coupling to gluons and to photons are affected as 
Sc s = ( ^c,A,(t.) + 2C 2 (v,)c i A,(t i ), 

fc, = ^,w + ^ f (,)-| t A(a (2.4) 

where 5q = Q — c^sm, C 2 (r) is the quadratic Casimir of the color representation, Tr(T a T b ) = 
C2(r)5 ab , and N(r) is the dimension (the number of colors) of the representation r. The 
functions Ai describe the 1-loop contributions of scalar, fermion, and vector particles to the 
triangle decay diagram. They are defined as 

Mr) = 4[/(r)-r], 

Mr) = 2^[(t-1)/(t)+t], 

Mr) = ^2 [3(2r-l)/(r) + 3r + 2r 2 ] , 

farcsinV^ r<l 
/ (r) = { x r ia , (2-5) 



T > 1 



1 

with Ti = m 2 h /Amf. Note that since f(r) ~ r + r 2 /3 for r < 1, one finds Ai(r) ~ 1 for the 
contribution of heavy particles (2m, ^> m^) . 



In the SM c g and c 7 arise from integrating out the top quark, and thus 

c 9 ,sm = A f (r t ) ~ 1.03 , c 7iSM = (2/9)c 9 , S m ^ 0.23 . (2.6) 
The invisible Higgs width in the SM is negligibly small, 

Cinv,SM ^ . (2.7) 

2.3 Partial Decay Widths and Branching Fractions 

With the help of the effective theory parameters, Cj, we can easily write down the partial 
Higgs decay widths relative to the SM values. Starting with the decays mediated by the 



lower-dimensional interactions in the first line of Eq. (2.1), we have 



bb 



n |2 F SM r ^ I |2 F SM r ^ I |2 F SM r |2 r SM (0 n\ 

c b\ 1 bb > 1 tt — \C T \ 1 TT , J- WW — \Cv\ 1 WWi 1 ZZ — \Cv\ 1 ZZi 



where, for m h = 125 GeV, the SM widths are given by Tf b M = 2.3 MeV, = 0.25 MeV, 



rffi w = 0.86 MeV and Tf§ = 0.1 MeV p]. Strictly speaking, Eq. gSb is valid at leading 



order. However, higher order diagrams which involve one q insertion leave these relations 



intact. Thus, Eq. (2.8) remains true when higher order QCD corrections are included. 
The decays to gluons and photons are slightly more complicated because, apart from the 
dimension-5 effective coupling proportional to c g , c 7 , they receive contribution from the loop 



of the particles present in Eq. (2.1). One has 
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-n _ \ C g\ -nSM -p _ I C 7 1 pSM icy Q \ 

l C 3,SM| l C 7,SM| 

where, keeping the leading 1-loop contribution in each case, one finds, 

c g = c g + c b A f (T b ) + c c A f (r c ), (2.10) 
c, = c,- 7 ^A v (r w ) + ^A f (r b ) + ^A f (r c ) + C ^A f (r T ). (2.11) 

Numerically, for m h ~ 125 GeV, A v (r w ) ^ 1.19 and A f (r b ) ~ -0.06 + 0.09z, so that 

c g ~ 1.03c g - 0.06c 6 , c 7 ~c 7 -1.04c v . (2.12) 



Consequently, c Si sm ~ 0.97, c 9ii sm ~ —0.81. The SM widths for that same mass are ~ 
0.34 MeV and ~ 0.008 MeV. 



In order to compute the branching fractions in a given channel we need to divide the 
corresponding partial width by the total width, 

Br(/i ii) = tix u = ^ . (2.13) 

1 tot 

The latter includes the sum of the width in the visible channels and the invisible width 



which, for m h = 125 GeV, is T inv ~ 1.2 x 10 3 |c in ^| 2 r|^. We can write it as 

r tot = |c tot | 2 r^, (2.14) 

where, for m h = 125 GeV, ~ 4.0 MeV, and 

ir I 2 r 

\r I 2 ~ \r, l 2 Rr SM -U \r,A 2 /Rr SM I Rr SM "l I 1 gl Rr SM I \r l 2 Rr SM I \r l 2 Rr SM I lnv 
\U to t\ — \Cb\ DT bb +\C V \ {DT WW + DT ZZ ) + 2 t>T gg + \C T \ tiT TT + \C C \ til cc + 

\Cg \ 1 tot 

ic I 2 r 

~ 0.58|c 6 | 2 + 0.2A\c v \ 2 + 0.09^^ + 0.06|cJ 2 + 0.03|c c | 2 + 

|as M | 2 

Typically, the total width is dominated by the decay to 6-quarks and T tot ~ c 2 , however 
this scaling may not be valid if the Higgs couples more weakly to bottoms (q, ^ 0.7), more 
strongly to gauge fields (cy > 1.4), or if it has a significant invisible width (cj nt) > 0.03). 

2.4 Production Cross Sections 

Much like the decay rates, one can express the relative cross sections for the Higgs production 
processes in terms of the parameters q. For the LHC and the Tevatron the currently relevant 
partonic processes are 

• Gluon fusion (ggF), gg — » ft+jets, 

• Vector boson fusion (VBF), qq — > hqq+]ets, 

• Vector boson associate production (VH), qq — > hV+jets 

The relative cross sections in these channels can be approximated at the leading order by, 
°" F - V ° 9? aVBF ~\c v \\ ^-kvl 2 . (2.15) 



^SM — is 12' n-SM ' v ' ' „SM 

a ggF l C 9,SM| a VBF a VH 



Using Eq. (2.15), we find that the total inclusive pp — » h cross section a to t, 



2^.SM l\% _ |2 j_ L |2^SM , \„ |2^SM 

(2.16) 



Is I2.-.SM /IS 12 i |„ |2„SM I \ r . \2„ 
&tot . . \ C g\ a ggF/\ C g,SM\ + \Cv\ °VBF+\ C V\ &VH 



°tot °ggF °VBF a VH 



is typically dominated by the gluon fusion process, and therefore it scales as a tot ~ c 2 . 



2.5 Event Rates at the LHC and Tevatron 

The event count in a given channel depends on the product of the Higgs branching fractions, 
the production cross section, and the selection efficiency. Experiments present the Higgs 
results as constraints on R (sometimes denoted as \i or p) defined as the event rates relative to 
the rate predicted by the SM. These rates can be easily expressed in terms of the parameters 



of the effective Lagrangian in Eq. (2.1) 



We first discuss the inclusive ATLAS and CMS searches in the 77, ZZ* and WW* final 
states. In those cases, the selection efficiency is similar for each Higgs production channels, 
to a good approximation. Thus, these searches constrain the inclusive rates R mcl defined as 
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(2.17) 



The approximations above hold assuming the Higgs production remains dominated by the 
gluon fusion subprocess. In our fits we use more precise expressions including the contribution 



of all leading production processes listed in Sec. 2.4 



Another category is associated production search channels targeting a particular Higgs 
production mode. One important example are the ATLAS and CMS studies in the 77 + 
2 jets final state where kinematic cuts on the jets were employed to enhance the VBF 
contribution [31 [7] . In the following we refer to these search channels as the dijet 77 category. 
Taking into account the selection efficiencies e, for different production channels, the dijet 
category searches constrain the rate R^ et defined as 



R5, 



e 99F\Cg\ 2 (Tgf F /\c g) SM\~ Z + e VB F\Cv\ Z 0-^BF + ^M^V^ Br(fe -)■ 77) 

e ggF°^ + tVBFO-y^p + e VH a VH Br SM (^ -» 77) 
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(2.18) 



where we used the approximated values given in Eq. (2.2). 



Finally, the searches in the bb final state at the Tevatron and the LHC target the Higgs 
produced in association with a (leptonically decaying) W or Z boson, therefore they constrain 



i?^ H defined as 



R 



VH 
bb 



a(pp — > Vh) Br(/i — > bb) 
osm(pp -> Vh) Br SM (h -> bb) 



CyCb 



a 



tot 



(2.19) 
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Figure 1: The combined signal strength (m = Ru and the corresponding error for the Higgs search 
channels used in our analysis. 



In summary, using Eqs. (2.8)-(2.19) we can express the observable event rates Ru in terms 



of the parameters of the effective theory denned by the Lagrangian (2.1). In the following 



we use the latest experimental determinations of Ru with the corresponding errors (assumed 
to be Gaussian) to identify the preferred regions of the parameter space. 



3 Combination and Global Fits 

Recently, the LHC updated the Higgs searches adding almost 6 fb _1 per experiment, collected 
at y/s = 8 TeV. Meanwhile, the Tevatron presented a new refined analysis of the full 10 fb _1 
data set. Here we focus on the following channels: the inclusive and the dijet tag h — > 77 
[31 E|, h ZZ* ->■ Al H U, h -)■ WW* -> 2l2v \9\ , and Vh -)■ Vbb M, E] channels. 
These are currently the most sensitive search channels for ~ 125 GeV. The combined 
central values for the signal strengths Ru in these channels and the corresponding errors are 
displayed in Fig. [T] In the following sections we use these bounds on Ru to constrain the 



parameters of the effective theory defined by (2.1). Throughout we assume = 125 GeV 
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and Gaussian statistics, and we do not take into account systematic effects which may be 
significant in the ZZ* channel and the dijet tag channel of the diphoton analysis. 

Since the ATLAS experiment did not provide the values of R^jj for = 125 GeV, we 
estimate the best fit values and error bars from the p- values in Fig. 11 of |2j. Assuming the 
probability distribution function follows Gaussian statistics, the best fit rate, /t, and error 
bar, a can found by inverting the observed and expected p- values (po)'- 

If a = dT^l - pg«P), ji/a = dT^l - p° bs ) (3.1) 

where 3>(x) = |(1 — Erf(x)) is the cumulative distribution function. Assuming the signifi- 
cance across each channel adds in quadrature, the rates from the 10 diphoton channels can 
be related to the rates of the 9 inclusive channels and the dijet channel. This can be used 



to calculate jljjjj and cr 77J j, after utilizing (3.1) to determine the other rates and error bars. 
For m h = 126.5 GeV, we find fi^jj = 2.9 ± 1.3, whereas ATLAS reports ft y7 jj = 2.7 ± 1.3, 
which we take as validation of our prescription for determining the rates using the reported 
p- values. 

With enough data from the LHC one could in principle perform a full seven-parameter 
fit, however for the time being we pursue a simpler approach. Here we assume c T = c& and 
Qnv = 0, and study the LHC and Tevatron constraints on the remaining parameter space. 
In this space, the best-fit points are 

\cy\ = 1.2l°;I, |cv| = 1.2l°;t \c g /c b \ = 0.818X \c b \ > 0.15. (3.2) 



Notice from (2.17)-(2.19) that x 2 mostly depends on the combination |c 9 /q,|, and thus does 
not strongly constrain c g and c& separately. The corresponding error bars are la, and have 
been derived assuming that the x 2 distribution around the minimum follows a Ax 2 distri- 
bution for 4 DOF. Relative to the SM point, A% 2 = Xsm ~ Xmm = 5-2, which is consistent 
with the global minimum at 74% CL. 

We also study the best-fit regions in new physics models where only two of the above 
parameters can be freely varied, while the remaining ones are fixed to the SM values. Sample 
results are displayed in Fig. [2] The results are shown for 4 different sets of assumptions about 
the Higgs couplings that can be realized in concrete models. In each plot we show the la 
constraints from the individual Ru while the "Combined" region corresponds to Ax 2 < 5.99, 
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Figure 2: The allowed parameter space of the effective theory given in Eq. (2.1), derived from the 
LHC and Tevatron constraints for nih = 125 GeV. We display the la allowed regions for the rates 
in Eqs. (2.17)-(2.19): i? 77 (purple), Rzz (blue), Rww (light grey), Rj-yjj (beige), and R b i (orange). 



The "Combined" region (green) shows the 95% CL preferred region arising from all channels. The 
crossing of the dashed lines is the SM point. The (^) corresponds to the best fit point. The top- 
left plot characterizes models in which loops containing beyond the SM fields contribute to the 
effective 5-dimensional hG^G^ and hA^ u A^ u operators, while leaving the lower-dimension Higgs 



couplings in Eq. (2.1) unchanged relative to the SM prediction. The red-line shows the trajectory 
5c.y = 2/95c g characteristic for top partners. The top right plot characterizes composite Higgs 
models. The bottom plots characterize top partner models where only scalars and fermions with 
the same charge and color as the top quark contribute to the effective 5-dimensional operators, 
which implies the relation <5c 7 = (2/9)5c g . In the bottom-right plot, the shaded region is 95% CL 
excluded by monojet searches at the LHC. 
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that is the 95% CL favored region. An interesting feature of all these plots is the presence 
of two disconnected best-fit regions [T^l [13] • 

The top left plot characterizes models in which loops containing beyond the SM fields con- 
tribute to the effective h G a ^ u G a ^ v and h A^A^ operators, while leaving the lower-dimension 



Higgs couplings in Eq. (2.1) unchanged relative to the SM prediction. Note that in this plot 
the band corresponding to i?^ H is absent; that is done for aesthetic reasons since -R^ H is very 
weakly independent of c g and c 7 . The data prefer negative contributions to c g (decreasing 
the ggF production rate) and to c 7 (increasing the diphoton decay width). An improvement 
of the fit with respect to the SM by A% 2 ~ 4 is possible for large 5c^/5c g , which would 
require a loop contribution from a particle with a large charge-to-color ratio [TH] {Qf ~ 2 
for the fundamental representation of SU(3) C ), or simply without color charge. Another 
possibility is a particle with Sc^/Scg ~ 0.25 (Qf ~ 1/2 for the fundamental representation 
of SU(3) C ) giving a very large negative contribution to the effective gluon coupling, tuned 
such that 5c g ~ — 2c 9i sm- Finally, a number of particles with different electric and/or color 
charges could be involved in such a way that their net contribution to 5c g (but not to 5c 7 ) 
approximately cancels. 

In the remaining 3 plots we fix <5c 7 = (2/9)5c g , which is the case in top partner models 
where only scalars and fermions with the same charge and color as the top quark contribute 
to these effective five-dimensional operators. In the top-right plot, the couplings to all the 
SM fermions, including the one to the top quark in the UV completed model is assumed 
to be rescaled by c&, producing the corresponding shift of c g and c 7 in our effective theory. 
Moreover, the coupling to W and Z is independently rescaled by cy. This is inspired by the 
composite Higgs scenario |20j . in which context a part of the parameter space with c& < 1 
and cy < 1 can be reached in specific models. The presence of two disconnected best-fit 
regions reflects the degeneracy of the relevant Higgs rates in the VV* and bb channels under 
the reflection q, — > — q,, which is broken only by _R 77 . A good fit is possible in the q, < 
region, although it may be difficult to construct a microscopic model where such a possibility 
is realized naturally. It is worth noting that the fermiophobic Higgs scenario, corresponding 
to q, = and cy — 1, is disfavored by the data (more generally, the fermiophobic line q, = 
is disfavored for any cy). 

The bottom-left plot demonstrates that the current data show a preference for a slightly 
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enhanced Higgs coupling to the electroweak gauge bosons, cy > 1. Several well-studied 
models such as the MSSM or the minimal composite Higgs (and more generally, models with 
only SU(2) singlets and doublets in the Higgs sector), predict cy < 1. If cy > 1 is confirmed 
by more data, it would point to a very specific and interesting direction for electroweak 
symmetry breaking [21J. 

In the bottom-right plot we relax one of our assumptions and allow for invisible Higgs 
decays. Opening the invisible channel reduces the visible rates, therefore it is in tension with 
the observations in the diphoton, VV*, and bb channels. Therefore, the data disfavor an 
invisibly decaying Higgs with Br inv > 20% unless there is a significantly negative 5c g (and, 
in the present case, the correlated negative contribution to <5c 7 enhancing the diphoton Higgs 
rate) in which case Brj ni , as big as 65% is allowed. Note that invisible Higgs production is 
directly constrained by monojet searches at the LHC [22J. As a result, interesting regions of 
the Br inv -5c g parameter space where the Higgs production cross section is enhanced compared 
to the SM are already excluded at 95% CL. 



4 Simplified Models 

In this section we discuss simple models where the interactions of the Higgs boson with 
matter may deviate from the SM predictions. In each case, we first map these models to 



our effective Lagrangian in (2.1), which facilitates an extraction of the observable rates. 
The models we consider introduce the minimal number of new degrees of freedom. We pay 
special attention to whether the new degrees of freedom allow for an enhanced diphoton 
rate, which is hinted by the data. We also study whether the models may improve the Higgs 
naturalness, i.e., whether they can cancel the quadratic divergent contributions to the Higgs 
mass induced by the SM particles. 

The study of such simplified models allows one to identify the required Weak-scale physics 
which can better account for the Higgs data. More complete models, which are typically 
constrained by additional experimental data can then be derived, and we postpone the 
detailed study of complete models to an upcoming publication. 
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4.1 Single Partner Models 

As a first exercise, consider a class of simplified models with only one new degree of freedom 
coupled only to the Higgs boson. The new degree of freedom, below referred to as the 
partner, could be a scalar S, a Dirac fermion /, or a vector boson p, carrying charge and/or 



color, and coupled to the Higgs as in Eq. (2.3) 



C = -c^htfS - cp-hff + c^hpfa . (4.1) 

V V V ^ 

Here appropriate index contractions are implicit for colored partners. For q = 1 the mass of 
the partner originates completely from electroweak symmetry breaking with a single Higgs, 
but we do not require this to be the case in general. For simplicity we assume in this 
subsection that the partner does not mix with the SM fields. This can be arranged, for 
example, by imposing a conserved Z 2 symmetry. We relax this assumption in the subsections 
below. 

Integrating out the partner affects the dimension-5 Higgs couplings to gluons and photons 
c g and c 7 , while keeping the remaining parameter in (2.1) at the SM values cy = q, = c c = 
c T = 1. The ratio c 9 /c 7 is determined by the electric charge and the color representation of 
the partner. As an illustration we present our results for the following three cases: 

• Scalar top partner. Color triplet, charge 2/3 scalar, contributing as 

5c g = C -^A s {ml/Am 2 s ) , (4.2) 
5c 7 = ^c s A s {m 2 h /Am 2 s ) . (4.3) 

The partner exactly cancels the quadratic divergence from the top quark for 



2m t ,A A\ 

c s = ^f- 4.4 

mi 



Fermionic top partner. Color triplet, charge 2/3 fermion, contributing as 

^9 = CfA f (m 2 h /Am 2 f ) , (4.5) 



2 



5c 7 = -CfAfiml/Amj) . (4.6) 



13 




Figure 3: Best fit regions in the Cj-mj plane, assuming = 125 GeV for the scalar top partner 
(top-left), fermionic top partner (top-right) and vector VF-partner (bottom- left). Shown are 
68% (darker green) and 95% CL (lighter green) regions. The dashed curves are for constant iCSr, 
Eq. (2.17), while the red curve is where a single partner is improving the naturalness of the SM, 



Eqs. (4.4), (4.7), (4.10). The bottom-right image shows the constraints for = 125 GeV, for 
top partner models, i.e. <5c 7 = 2/95c g . The three bands show the la allowed regions for RY. , 
, and R^z channels. The three curves show the theoretical predictions as a function of 5c g for 
each channel. Only 3 channels are shown, but all channels are included. The green shaded region 
shows the 95% CL experimental preferred region. 
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The partner exactly cancels the quadratic divergence from the top quark for 

2m t ,A ^ 

Cf = 7 / \ , (4-7) 

rrif ( rrif + ./2m ( 2 + mj J 

and, in that case, in the limit m/ 3> m t one has 8c g = 9<5c 7 /2 rs — m|/m,y. 

W prime. Color singlet, positively charged massive vector, contributing as 

5c, = , (4.8) 
7 

<$Gr = --c p 4,(m£/4m?) . (4.9) 



The partner exactly cancels the quadratic divergence from the W boson for 



rn 



w 



(4.10) 



The best fit 68% CL and 95% CL regions for these 3 examples are shown in Fig. |3j In 
each case, the model is defined by 2 parameters: the coupling q and the mass rrii, but only 
one combination influences the effective theory parameters. Consequently, the best fit region 
in the Cj-m.j plane corresponds to a line of minimum A% 2 . For a scalar, the best fit occurs 
for a large negative coupling c s , which is incompatible with naturalness in this simple set-up. 
To improve the fit to the Higgs data with scalar partners one needs at least two of them 
(as in supersymmetric theories) and with a large mixing, so that the sign of coupling of the 
lighter partner to the Higgs is flipped due to the mixing angle. Note that in that case the 
lighter scalar actually worsens the quadratic divergence of the SM, and the restoration of 
naturalness is postponed to the higher scale where the heavier scalar intervenes. 

For the fermionic top partner case the situation is different as a negative coupling is 
consistent with naturalness. Therefore a single fermionic top partner may improve the fit 
to the Higgs data and the naturalness at the same time provided the partner is light, with 
mass in the narrow range of 95 GeV < m/ < 115 GeV. Such low masses are likely to be 
excluded by direct searches and precision measurements in specific more complete models. 
In the vector W partner case, there are 2 separate best fit regions with the same Xmm : one 
where c p is relatively small and positive, where the partner interferes constructively with the 
SM W boson, and the other where c p is large and negative, so that the partner "overshoots" 
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Figure 4: Constraints for the Simplest Higgs model (left) and the Twin Higgs model (right) 
assuming = 125 GeV. The three bands show the la allowed regions for R^ b H , R™ 1 , and R™z' 
channels. The three curves show the theoretical predictions as a function of £ for each channel. 
Only 3 channels are shown, but all channels are included. The green vertical lines show the 95% 
CL experimental preferred region. 



the SM W contribution. Only the latter region can overlap with the curve where the partner 
exactly cancels the quadratic divergence from the W boson. 

To summarize, for all models an enhanced diphoton decay rate can be obtained. In the 
scalar and fermionic case one needs a sizable negative coupling q, while in the vector case a 
moderate positive coupling suffices. The diphoton rate is indicated by the constant contours 
shown in Fig. |3} For the top partner case, the fit to data is improved relative to the SM as 
A X 2 = Xmin - Xlu = 3-8 for only 1 dof. 



4.2 Composite Higgs Models 

The single fermion partner model described in the previous subsection is a special case of 
models with a more general set of couplings of fermions to the Higgs. In particular, consider 
a vector-like pair of top partners, T and T c , which interact with a pseudo-Goldstone Higgs. 
There are several possibilities for the interactions, depending on the type of model under 
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considerations. Here we consider two possibilities: 

- C top = yf sm(\H\/f)tX c + yf cos(|tf|/ f)Tt c + M'TT C + h.c. . (4.11) 

where either X c = t c or X c = T c . The former occurs, e.g., in the Simplest Little Higgs model 
with an [SU (3) / SU (2)] 2 coset [23] or in the minimal composite Higgs with SO(5)/SO(4) 
coset structure [25J. The latter case, is encountered in the left-right Twin Higgs model [24J. 

For sufficiently large M', the heavy partners can be integrated out and the low energy 
couplings to gluons and photons are found to be, 



Cg Cry 



i - iii x c = t c 
, 1 t > ( 4 - 12 ) 

where £ = v 2 /f 2 . Similarly, for the models above one finds, 



c v = c 6 = Vl-C/2, (4.13) 

Thus we see that the couplings of the pseudo-Goldstone Higgs may exhibit a universal or 
non-universal suppression which depends on the single parameter, £, and are independent 
on the specific details within the top sector (masses and mixing, for example). 

In Fig. |4j we show the production and decay rates of the above models as a function of 
£. We find that for a 125 GeV Higgs boson, £ is constrained at 95% CL to be £ < 0.4 in the 
universal suppression case, and £ < 0.2 in the non-universal case. 

4.3 Dilaton 

Here we study the hypothesis that the 125 GeV resonance discovered at the LHC is a dilaton 
[261 [14"] . that is to say a (pseudo)-Goldstone boson of the spontaneously broken conformal 
symmetry. That implies that it couples to the trace of the energy-momentum tensor of the 
SM matter, 

C = c£ Um 2 w W+W- + m 2 z Z,Z, - ^ m f ff ] ~ 0T£ (4.14) 
V V feSM J 

In other words, the lower-dimension couplings are the same as for the SM Higgs bosons, up 



to an arbitrary overall rescaling factor c^. Mapping to the effective Lagrangian (2.1), we 
have, 

c v = c b = c T = c c = c^. (4.15) 
17 




Figure 5: Left: The difference between the x 2 °f the dilaton model and the x 2 °f the SM as a 
function of the parameter ca. At the best fit point around q, ~ 0.27, the x 2 °f the dilaton model 
is larger by 5.2 units compared with the SM, indicating that the dilaton model always fits the data 
worse than the SM. Right: The favored region at 68% CL (Darker) and 95% CL (Lighter) for the 
125 GeV resonance being a mixture of the SM Higgs and a dilaton. The best fit occurs along the 
line a = corresponding to a pure SM Higgs. The dashed lines show contours of constant R~~. 
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The dimension-5 coupling do not however follow the overall rescaling because they are af- 
fected by integrating out new heavy degrees of freedom needed to restore conformal invari- 
ance at high energy Using conformal symmetry, these UV contributions to c g and c 7 can be 
related to the QCD and electromagnetic beta functions: c g = — c^36q CD /2, c 7 = — c^6q M /8. 
Since above the scale of the top mass one has 6q CD = —7 and 6q M = H/3, integrating out 
the top quark one finds 

21 11 

6c 9 = y c </> + ( c <t> - !) A A T t) , 5c 7 = ~^ c <t> + ~ x ) A /( r *) • ( 4 - 16 ) 

Thus, for > 0, the effect of the dilaton is to increase the H — > 77 width, which is favored 
by the data, but at the price of increasing the ggF production rate, which is disfavored. We 
therefore find that at the best fit point at ~ 0.27, the x 2 °f the dilaton model is larger 
by 5.2 units than Xsmi as can be seen in the left panel of Fig. [HJ Thus, the simple dilaton 
interpretation of the 125 GeV resonance is not favored by the data. Allowing for negative 
does not change that conclusion, as the rates are symmetric under — > — c^. 
More generally, one may consider the dilaton mixed with the SM Higgs boson, with the 
mixing angle denoted by a. In that case the effective theory parameters are given by 

21 

5c g = — sin a + (c^ sin a + cos a — 1) Af(r t ) , 

5c 7 = — — c<^ sin a + (c^ sin a + cos a — 1) Af{r t ) . (4-17) 

The best fit regions in the a-c^ plane are shown in the right panel of Fig. [5j Again, the best 
fit is obtained along the SM line a = 0. 

We conclude, therefore, that at the moment, there are no hints of a dilaton nature in 
the 125 GeV resonance, although a large Higgs-dilaton mixing angle cannot be excluded at 
present. 

5 Extended Higgs Sectors 

We next turn our attention to extended Higgs sectors with one or more additional scalars 
beyond the 125 GeV Higgs. As we show below, several scenarios are possible, allowing for 
an enhanced diphoton or ZZ and WW rate. 
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Figure 6: Favored region at 95% CL, of the 2HDM in the a — tan/3 plane. For sin a < the 
favored region is concentrated around the decoupling limit (dashed line), a = (3 — ir/2 where all 
couplings are SM-like. Conversely, for sin a > the favored region lies around a = (3 + ir/2 where 
the top Yukawa coupling is SM-like, but cy = —1. 

5.1 The Two-Higgs Doublet Model 

Consider the Type II 2-Higgs Doublet Model (2HDM), with 2 Higgs doublet fields H u , H d : 



(5.1) 



s/2 J ' " a V y/2 ' d 
The doublet H u couples to up-type quarks, and H d to down-type quarks and leptons. The 
ratio of the two VEVs is tan /3 = v u /v^ and by convention < /3 < ir/2. The CP-even mass 
eigenstates are mixtures of the neutral components 







[hi 


H 



cos a sina j I Re(H^) 
— sin a cos a I \Re(i?°) 



(5.2) 



with a the mixing angle. We identify h as the 125 GeV Higgs. The tree level Higgs couplings 
to the fermions and vectors are given by 



sin a cos a 
cy = sm(p — a) , Cb = — - — ^ , c g = 9/2c ~ 



cos (3 



sin/3 



(5.3) 



where the effect on c g arises because of the modified Higgs coupling to the top quark. 
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Figure 7: Left: The favored region at 68% CL (Darker) and 95% CL (Lighter) for = 125 GeV 
in the two scalar model, with m ( - 2 3> , and assuming the decoupling limit of the 2HDM. The 
dashed lines show contours of constant i? 77 . Right: Allowed regions at 95% CL obtained for the 
2HDM + two stops model. The different colors correspond to different values of tan /3. The regions 
are: tan/3 = 5 (yellow), tan/3 = 10 (green), tan/3 = 25 (blue) and tan/3 = 50 (purple). The best 
fit points corresponds approximately to the decoupling limit. 

In Fig. [6] we show the constraints in the a-tan plane. In accordance with current direct 
bounds [27], we assume that the charged Higgs is heavy enough so as to contribute negligibly 
to c 7 , and tan /3 > 0.3, so that the top Yukawa coupling does not run to a Landau pole at 
fi ~ TeV. The best fit approximately corresponds to the decoupling limit a = (3 — n/2 
(tana = — cot /3) where all couplings are SM-like, in particular cy ~ 1. The minimum x 2 is 
roughly the same as in the SM, A% 2 = xLin ~ Xsm — 0- Another favored region is for a > 0, 
where cy is still close to 1 and the sign of q, is nipped. 

5.2 Simplified MSSM 

The next example we consider is a simplified model of the MSSM: two Higgs doublet plus 2 
stops defined as scalars with the same color and charge as the top quark. The Higgs doublets 
are defined as in the previous section, but now a G (— 7r/2,0). Consider the stops t,t c with 
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the mass terms of the form 



- C stop = \t\ 2 (m 2 + y 2 t \H u \ 2 ) + \f\ 2 (fh 2 c + y 2 t \H u \ 2 ) + y t \H u \X t (tt c + h. 



•c) , (5.4) 



where y is the top Yukawa coupling. This is equivalent to the stop sector of the MSSM 
neglecting the (sub-leading) D-terms contribution to the stop masses. The left-handed and 
right-handed stops mix in the presence of X t , which in the MSSM is given by X t = \A t — 
\x cot j3 1 . 

We begin by considering the decoupling limit, ^> m^. In that case, the change in 
rates are controlled by the stop spectrum. Denoting the two mass eigenvalues by m^., and 
the left-right mixing angle by 8t, one has 

m t X t = i (m| - m|) sin 20*, (5.5) 

where, by convention, < m t - 2 . For 3> m^/2, integrating out the stops shifts the 
effective dimension-5 operators as 



c„ Ov If m+ m 2 m 2 X? , 

c 9 ,sm c 7iSM 4 \ml mj 2 mfmf 2 



For zero mixing, the stops always interfere constructively with the top contribution (and de- 
structively with the IV-contribution) to c 7 . Assuming that both stops are heavier then about 



100 GeV, we see from Eq. (5.6) that the contribution to c g is bounded above by Sc g < 1.5. 
For large X t , the sign of the contribution from stops can flip, and a significant enhancement 
of the diphoton width is possible. Since X t needs to satisfy the bound 2m t \X t \ < (m~ —m~), 
it can be large provided the two eigenstates are split. 

In Fig. [7] we illustrate the impact of the LHC Higgs data on the parameter space of 
the simplified MSSM model. The left plot shows the preferred region in the m^-Xt/m^ 
plane, assuming the decoupling limit and m t ~ 2 3> so that the heavier stop eigenstate 
does not contribute to the effective operators. In this case, the Higgs data only constrains 



one combination of parameters, that is c g in Eq. (5.6). As can be seen in Fig. |7| , the data 
strongly prefers large mixing, where the mixing dominates the stop contribution to c g and 
c 7 . The best fit contour corresponds to 5c g = —1.78 and A% 2 = xLm ~ Xsm = 3-8 for 1 
DOF, thus the data strongly favors the top partner scenario. On the other hand, the no 
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Figure 8: The best fit regions in the a-5 7 plane, assuming 777^=125 GeV and the singlet coupled to 
top-like (left) and colorless (right) vector-like fermions. Shown are 68% (darker green) and 95% CL 
(lighter green) regions. The cross corresponds to the best fit point. The dashed lines are contours 
of constant i? 77 . 

mixing scenario, X t = 0, is strongly disfavored, since 5c g > 0. Indeed in that case, the stop 
contribution can only enhance the gluon fusion cross section while decreasing the diphoton 
branching ratio, which does not fit well the LHC data. 

Going back to the more general situation where the Higgs sector is away from the decou- 
pling limit, the right plot of Fig. [7]shows the constraints from the Higgs data in the Sc g — sin a 
plane for different values of tan /3. Although the parametrization is motivated by the MSSM 
with light stops, it is also applicable to any 2HDM model with top partners. There are 3 
independent parameters which are constrained by the Higgs measurements: tan/3, sin a, and 
5c g . Even with the additional degrees of freedom, the Higgs fit to the data is not improved 
relative to taking the decoupling limit. This is simply a consequence of the fact that the 
data strongly points to cy ~ 1, which can only be achieved in the decoupling limit. 
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5.3 The Doublet- Singlet Model (Enhanced c 7 ) 

The current CMS and ATLAS data point to a somewhat enhanced diphoton rate. One simple 
scenario which may allow for such an enhancement is the case where the Higgs doublet mixes 
with a singlet, thereby affecting the lower-dimension parameters of the effective Lagrangian 
(2.1). More precisely, let us assume that the SM Higgs doublet H = {0,v + h )/V2 mixes 



with a real scalar, <p, which couples to new vector-like charged fermions ipi, 

C D -\m 2 y - Ktp\H\ 2 - V Mi(l + -<p)AA • (5.7) 

i 

The two mass-eigenstate Higgs scalars, h and H°, mix with one another with an angle a 
defined by 

h = ho cos a + if sin a , H° = — h sin a + (p cos a , (5.8) 

and we assume h to be the Higgs at ~ 125 GeV. The mixing angle suppresses the 
couplings of h to WW, ZZ and SM fermions by cos a relative to the SM. On the other 
hand, h receives additional, possibly positive, contributions to the couplings to F^F^ and 
G a ^ v G a ^ v from loops of the fermions ipi, resulting in 

c g — c g,SM cos a + 5c g sin a , (5.9) 
c 7 = c 7i sm cos a + 5c.y sin a . (5.10) 

Here the "singlet" one-loop couplings 5c g and 5c 7 are given by 

5c g = J2 2C 2(n)kA f (n) , 5c 7 = \NiQl\iA f {Ti) . (5.11) 

i i 

The fermion representations and couplings, A«, can always be chosen to obtain any desired 
correction to 5c g and 5c 7 . 

Mapping to the effective theory one finds, 

Cy = Q, = c c = c T = cos a , (5-12) 

c g = cos a A f{r t ) + 5c g sin a , (5.13) 
2 

c 7 = - cos aAf(T t ) + 5c 7 sin a . (5-14) 
9 

Accordingly, the double-singlet model has three free parameters: a, 5c g , and 5c 7 . Two 
examples of motivated ansatzes are £c 7 = (2/9) 5c g (when ipi have quantum numbers of the 
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top quark), and 5c g = (when ipi are color singlets). The fits in the a-5c 7 plane for these 
two ansatzes are shown in Fig. S 1 Only <5c 7 > is shown as the region with <5c 7 < is 
equivalent upon a — > —a. The best fit is obtained for a negative mixing angle, which allows 
for an enhanced Higgs diphoton rate, (as demonstrated by the constant contour lines). On 
the other hand, the mixing angle is required to be sufficiently small in order to avoid strong 
suppression of the WW and ZZ rates. 

Finally, we note that the doublet-singlet model predicts an additional resonance which 
may look at colliders much like the Higgs field, but with suppressed couplings. While its 
mass remains a free parameter, the existence of such a state may therefore allow one to 
confirm or exclude this model in the near future. We postpone further details of this model 
for future work. 



5.4 The Doublet-Triplet Model (Enhanced c v ) 

If the Higgs sector contains triplets or higher representations under SU(2)w then the Higgs 
coupling to H^and Z bosons can be enhanced. The Georgi-Machacek (GM) model [281 [29]. 130] 
is one example that contains Higgs triplets, is renormalizable, and does not introduce large 
violations of the custodial symmetry. 

The Higgs sector of the GM model contains the usual Higgs doublet H transforming as 
2i/2 under SU{2)w x U(l)y, a real triplet (f> transforming as 3o, and a complex triplet A 
transforming as 3i. This field content forms irreducible representations of the approximate 
global symmetry group G = SU (2) L x SU(2) R : the doublet can be put in $ forming the (2, 2) 
representation of G, while the triplets can be put in \ forming the (3, 3) of G. The vev v = 246 
GeV is distributed between the doublet and triplets as (<&) = v 0/3/2x2, (x) — vs /3hx3- Since 
the T parameter is protected by a built-in custodial symmetry, the triplet vevs parameterized 
by sp = sin (3 can be 0(1) without conflicting phenomenology. The vevs of $ and x break 
G down to SU(2)v referred to as custodial isospin, under which $ decomposes as singlet 
H(2) and triplet GW), while x decomposes as singlet H@), triplet G(3), and quintuplet Q. 
Of the custodial triplets there are only three physical states, A^jA , while the Goldstone 
bosons, G, get eaten by W and Z. Here A and G are defined by, G^ = cpG a — spA a , 

2 We thank Matthew McCullough for pointing out an error in the earlier version of this plot. 
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G/U = spG a + cgA a , where a = ±, 0. The isospin singlets also mix, 

H{2) = c a h - s a H , = s a h + c a H , (5.15) 

where the angle a depends on the details of the Higgs potential. We identify the 125 GeV 
resonance with h. Further technical details about the model are postponed until Appendix |A| 
The GM Higgs potential in the custodial limit contains 7 parameters: 2 masses and 5 



quartic couplings, see (A. 9). They can be related to 7 observables: the known VEV v, the 
2 mixing angles a, (3, and the masses of the 2 isospin singlets, = 125 GeV and m#, the 
isospin triplet, m^, and the isospin quintuplet, uiq. In our convention the Higgs is SM like 
for a — [3 — 0. The light Higgs phenomenology is thus determined by the 4 parameters a, 
(3, rnq and m^, the former two affecting the lower dimensional Higgs couplings, the latter 
affecting the dimension-5 couplings to gluons and photons. Mapping to the effective theory 



in (2.1 ) we find, 



c v = c a c p + ^/8/3s a sp, c b = c c = c T = —, (5.16) 

^ 9 - c b A f (r t ), c 7 = 2 -c b A f (r t ) + 9 -^A s (t a ) + ^ s (r Q ). (5.17) 

Note that the coupling to W and Z can be enhanced over the SM value up to the maximal 
value of cy = a/8/3, which is a distinct feature of models that contain custodial quintu- 
plets under the custodial symmetry. The coupling to the fermions can also be enhanced or 
suppressed compared to the SM Higgs, depending on a and /3. The c g and c 7 differ from 
the SM value because of the modified Higgs-top coupling while c 7 also receives additional 
contributions from integrating out the charged Higgses Q ++ ,Q + and A + . The couplings 



9hA*A and ghQ*Q are given in (A. 12) and (A. 13). 



The best fit regions in the s a -sp plane fixing rug = = 300 GeV are shown in the left 
panel of Fig. |9j The Higgs data are consistent with s@ significantly larger than zero, that 
is with sizable triplet vevs. The best fit regions are concentrated along \s a \ ~ sp where the 
Higgs coupling to fermion is near the SM value, but the coupling to W and Z is enhanced, 
leading to an enhanced diphoton rate. Note that the best fit regions are well within ranges 
allowed by the other constraints (see also Ref. (31] for a recent phenomenological study of 
the GM model). The constraints on tua as a function os coming from b — > 57, B s ^ — B s ^ 
and t — > bA + decays are shown in the right panel of Fig. |9j We have also checked that 
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Figure 9: Left: Best fit regions at 68% (darker green) and 95% CL (lighter green) in the 
(sin a, sin /3) plane, setting nih = 125 GeV and toq = ttia = 300 GeV. The SM point (cross) 
compared to the \ 2 minimum (red dot) corresponds to x|m ~~ Xmm = 4-4. The dashed lines show 
contours of constant Rjj. Right: The 95% CL allowed regions for sin/3 and obeying con- 
straints from B s ^ — B St d mixing (blue band with solid contour), from 6—7-57 ( re d band with 
dashed contour) and from t — > bA + decays (small regions inside dotted curves are excluded). 

the ttiq = 300 GeV is not constrained by doubly charged Higgs searches at CMS [32J, 
ATLAS [33J, CDF [31] and DO [35], and from searches for anomalously large production of 
multi-lepton final states [35] (similar sensitivity is expected from SUSY searches that require 
same sign leptons with missing jet and MET [371 EHl 139]). 

6 Conclusions 

The recent discovery of a Higgs boson at 125 GeV has important consequence for any new 
physics scenario. In this paper, we studied those implications using the Higgs searches 
reported by ATLAS, CMS and the Tevatron in currently the most sensitive channels. We 
derived the constraints on the parameters of the effective Lagrangian describing, in a very 
general fashion, the leading order interactions of the Higgs particle with matter. Overall, the 
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point in the parameter space corresponding to the SM Higgs boson is well consistent, within 
74% CL, with the Higgs data. The data show a preference for models leading to an increased 
rate in the diphoton search channels, but at this point this is not statistically significant. 
We note that the current data are still consistent with a sizable invisible branching fraction 
of the Higgs: Bi inv < 20% if the Higgs production rate is the one predicted by the SM, and 
Br im , < 65%, if new large contributions to the effective Higgs coupling to gluons and photons 
are present. 

As a second step, we mapped a number of new physics models affecting Higgs phe- 
nomenology on to the effective Lagrangian, placing constraints on the relevant couplings. In 
particular, we study the simplified version of the MSSM, composite Higgs models, dilaton, 
a two-Higgs doublet model, a doublet-singlet model and a doublet-triplet model. In each 
case we identified the region of the parameter space that can improve the quality of the 
fit by enhancing the rate in the diphoton channel. If that feature persists with more data, 
indicating a departure from the Standard Model, new physics models would require special 
structure and additional charged particles, likely within collider reach, in order to explain 
the data. 

Note added: While this work was being completed the analyses [ID] appeared that 
overlap to some extent with our results. 
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A Further details on the Georgi-Machacek model 



We first introduce some notation. The J = 1/2 SU(2) generators are T a = a a /2, where a a 
are the usual Pauli matrices, while T" are the SU(2) J = 1 generators in charge eigenstate 
basis, 

/ i \ 

-i -i 
\ % / 



7? = 



1 



/ -1 o\ 



v/2 



-1 1 
\0 10/ 

We also defined the matrix P 3 



Tl = 



V2 



(l o \ 



Tf = 







(A.l) 



/ 1 \ 



-Fa = 



(A.2) 



1 

Vi o o/ 

that satisfies P 3 (T a )* = — T a P 3 . Thus, if ip transforms as a triplet, then ip = P3^* transforms 
as a triplet too. 

In the GM model, the electroweak doublet H and the triplets 0, A can be collected into 
(2,2) and (3,3) representation under the global custodial G = SU(2) L x SU(2) R : 



*=(hh), x = (A <\> a) 



(A.3) 



Here H = ia 2 H* and A = P3A*. In this notation the {/(l)y is realized as the T 3 generator of 
SU(2)r. The vevs of $ and x break G down to the diagonal SU (2) referred to as the custodial 
isospin, under which $ decomposes as singlet if (2) and triplet (7(2), while x decomposes as 
singlet if ( 3 ), triplet (2(3), and quintuplet These fields are embedded into the doublet and 
triplets as, 



H 



iG+ 2) 



( 







\ 



2V2 1 V 3^0) V 3 



1- \l *H( 3 ) - \ I 



±(v Cf3 + H {2) -iG {2) ) 

( Q ++ 



\ 



75(^ + ^3)) 



(A.4) 



(A.5) 



The kinetic terms in the Higgs sector are given by 

C kin = [D^D^] + ^Tr [D^D.x] , 



(A.6) 
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with the covariant derivatives 

r 3 

2 ' ~ a 2 
D, X = d, X ~ igTtKx + ig'B^xTf, (A.8 



= - ig°-L a ^ + ig'B^*-, (A.7) 



One can check that Eq. (A. 6) leads to the kinetic terms for the scalars that are diagonal 
and canonically normalized, and it also leads to the usual SM gauge boson masses with W 
and Z kinetically mixing with Goldstone bosons G and G°. The most general custodially 
invariant and renormalizable potential is [29] 

7T1 2 VP? 9 

V =^Tr [$t$] + _* T r [ X ] X ] + Ai (Tr [$t$]) + A 2 Tr Trf^x] 
+ A 3 Tr[ X t xx t x] + x 4 (Ti[ X j x}) 2 ~ A 5 Tr Tr^T^]. 
The quark Yukawa terms only involve the Higgs electroweak doublet, 

£ = — % + £T( 2 )) QlQr -> -m, ( 1 + — - - — — J g L g H , (A.IO) 

where m g = y q vcp/\/2. In the lepton sector there the Yukawa couplings to both the Higgs 
doublet H and the triplet A are possible. The coupling of the doublet is of the same form as 
for the quarks, C = —yeHLe c +h.c, while the coupling to the triplet, C = —\iA a La a L+h.c., 
also contributes to Majorana neutrino masses as soon as sp ^ 0. Given the smallness of the 
neutrino masses the Yukawa coupling constants are exceedingly small for realistic Sp, and 
can be neglected for our purposes. 

The h — > 2j decay is affected by the charged component of the custodial isospin triplets 
and quadruplets running in the loop. Keeping only the trilinear couplings of h to the triplet 



and quintuplet in the potential (A. 9) one gets 

C = -g hA *Am 2 A h (2A + A- + A 2 ) - g hQ * Q m 2 Q h (2Q++Q" + 2Q + Q~ + Q 2 ) , (A.ll) 
where 

ghA*A = (c a C/3 + y/8/3s a s p ) + (m 2 h /m 2 A ) (2y/^s a + 3c a s 3 p)/ (6c p s fi ) , (A. 12) 

g hQ * Q = a/273 (2 + m 2 h /m 2 Q )s a /sp + (m 2 A /m 2 Q ) cp(-2y/6cps a + 3c a sp)/sp. (A. 13) 
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Let us now review the constraints on masses of the extra scalars, H, A and Q. The heavy 
Higgs is constrained by direct searches at the LHC for the SM Higgs decaying to WW and 
ZZ. The bounds are avoided for m# > 600 GeV, while for smaller values sp is constrained 
[3T] . The charged A + boson couples to fermions through the Yukawas which in the mass 
eigenstate basis are 

Ant = i{u L V C KMYdd.R)SfiA Jr + i(u R Y u V ( l KM d L )spA + + h.c, (A.14) 

where Y u ^ = di&g(m Ut d)\/2/v and Vckm is the CKM matrix. The FCNCs are generated at 
1-loop. The most constraining are b — > 57 and B^ s — B^ s mixing with A + and top in the 
loop. To obtain the bounds we use LO matching onto effective weak Hamiltonian at \xy/. For 
b — > 57 we use the LO equations (53-56) in [U] with the replacement A u ^ = ±sp/ {v\J2Gf) 
(the NLO results are also available [HJ H2]), the parametrization in Eqs. (46-49) of |43] for 
the effect on 6 -)■ 57, the experimental value of BriB — > X^) = (3.55 ± 0.26) ■ 10~ 4 [H] 
and the SM prediction Br(B X s7 )sm = (3.15 ± 0.23) • 10" 4 [45] (for the estimate of the 
irreducible error, see [16]). The resulting bounds on sp and are shown in Fig. [9] as a red 
band. For constraints from B^ s — B^ s mixing we used the LO matching expression in Eq. 
(54-56) of [47] with the replacement 1/ tan/3 — > s@. The ratio of NP to SM contributions 
hd yS = ^i2 P /^i2 M = (2Swh + Shh) I Sww is the same for Bj — Bj and B s — B s systems since 
the Georgi-Machacek model is an example of an MFV extension of the SM. We therefore 
use the results of the fit where MFV constraints are imposed on the mixing amplitudes, 
hd, s = _ 0.08^q;q3 9 [IS] (this is Scenario II in the terminology of [48J) leading to constraints 
on sp and itla shown in Fig. [9] as a blue band. Note that in Georgi-Machacek model hd )S is 
always positive. Since the central value of hd, s is negative experimentally (but well within 
one sigma range) this makes the constraints slightly more stringent. 

Finally, the direct searches relevant for Q constraints are direct searches for H ++ , H + — > 
rv and SUSY searches with leptons and MET. In the doubly charged Higgs searches from 
CMS [32], ATLAS [33], CDF [31] and DO [35] the dominant decay modes were assumed to 
be Q ++ — > £ + £ + and Q + — > £ + v. In Georgi-Machacek models the dominant decays are not 
into leptons but into gauge bosons, which make the searches not sensitive to Q at present. 
Similarly we have checked that the search for anomalously large production of multi-lepton 
final states is not sensitive to mq ~ 100 GeV [36] and similar (in) sensitivity is expected from 
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other SUSY searches [37] EU [39]. 

More constraining are the charged Higgs searches in the H + — > tv channel, both in CMS 
and ATLAS [SO]- These are only relevant, if A + is lighter than t. In this case the 
branching ratio for Br(t — > bA + )Br(A + — > r + u) is bounded to be below between 4% for 
rriA = 90 GeV and 1% for rriA = 160 GeV at 95 C.L. [50] (the limits are comparable for CMS). 
The branching ratio Br(t — > A + b) ~ s 2 ^ in Georgi-Machacek model, if we neglect the effect of 
phase space for clarity. On the other hand Br(A + — > r + v T ) : Br(A + — > cs) : Br(A + — > cb) ~ 
m 2 . : m 2 c : \V c b\ 2 m 2 = 0.06 : 0.03 : 0.91, with masses evaluated at the scale fi ~ rriz- This 
means that A + — > t + u t is a subdominant decay mode with only ~ 6% branching ratio. The 
ATLAS search thus excludes at 95%CL a small region of (m^, s@) parameter space around 
rriA ~ 120 GeV, see Fig. [9] (CMS is just about to become constraining). 
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